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d ' Abstract 

CJ . The motivation for the treatment of intrabeam scattering theory 

*55 ! given in this paper was to find results which would be convenient 

J^"* for computing the intrabeam scattering growth rates for particle dis- 

O ■ tributions which are more complicated than a simple gaussian. It 

was shown by A. Piwinski that beam growth rates due to intrabeam 
scattering can be expressed as a multidimensional integral [1]. It was 
!> ■ pointed out by J. Bjorken and S. Mtingwa [2] that the reduction of the 

multidimensional integral to a 3-dimensional untegral is made easier 
| by writing the integral so that its relativistic transformation proper- 

•■O ■ ties are more obvious. The starting point in [2] was a result from the 

treatment of the two body scattering problem in relativistic quantum 
£^ \ theory . In this paper the starting point is the relativistic transfor- 

""^5 ' mation properties of the scattering cross section which may be a more 

familiar starting point. The resulting expression for the multidimen- 
sional integral is simpler to reduce. In addition, the results do not 
depend on the particular form of the Coulomb cross section that was 
Qh used in [2] and are valid for any collision cross section. 

> '■ 

^ ! 1 Introduction 

a, 

The motivation for the treatment of intrabeam scattering theory given in 
this paper was to find results which would be convenient for computing the 
intrabeam scattering growth rates for particle distributions which are more 
complicated than a simple gaussian. It was shown by A. Piwinski that beam 
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growth rates due to intrabeam scattering can be expressed as a multidimen- 
sional integral [1]. It was pointed out by J. Bjorken and S. Mtingwa [2] that 
the reduction of the multidimensional integral to a 3-dimensional untegral 
is made easier by writing the integral so that its relativistic transformation 
properties are more obvious. The starting point in [2] was a result from the 
treatment of the two body scattering problem in relativistic quantum theory 
. In this paper the starting point is the relativistic transformation properties 
of the scattering cross section which may be a more familiar starting point. 
The resulting expression for the multidimensional integral is simpler to re- 
duce. In addition, the results do not depend on the particular form of the 
Coulomb cross section that was used in [2] and are valid for any collision 
cross section. The final result is given by Eq.(14), which can be used for com- 
puting the intrabeam scattering growth rates for particle distributions which 
are more complicated than a simple gaussian. 

2 Transformation properties of the cross sec- 
tion 

The cross section, a, which describes the the scattering of particles with the 
momentum p 1 from the target particles with momentum p 2 is first defined 
in the CS (coordinate system) where the target particles are at rest, p 2 = 0. 
In a scattering event, the particle momenta change from pi,p 2 to p'i,p 2 - As 
we are assuming that both momentum and energy are conserved , the final 
momenta, p[,p 2 are determined by the direction of p[ which is indicated 
by the unit vector p[. In this CS where p 2 = 0, a is defined so that the 
number of incident particles which are scattered by the target particles with 
momemtum p 2 which are in the volume element, d 3 x, in the time interval dt, 
into the solid angle dQ! corresponding to the direction p[ is given by 

SN = adQ'p 1 (x)v 1 p 2 (x)d 3 xdt (1) 

where pi(x),p 2 (x) are the density functions and vi is the velocity of the 
incident particle. 

Now let us go to a CS where p 2 ^ 0. In this CS, a is defined by requiring 
adQ' to be invariant, that is to have the same value in all coordinate systems. 
A simple way to find the relationship between SN and a in this CS is to write 
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SN as (see [3]) 



5N = ad n>'^tl^±F( Pl ,p 2 )d 3 xdt (2) 
7i 72 

where p\{x)/^i is an invariant as it is just the density function for particle 1 
in the CS where p\ = 0. Similarly for p 2 {x)/ r y 2 . If one can find an invariant 
F(pi,p 2 ) which for p 2 = gives F = 7ii>i, then this expression for SN gives 
the correct result when p 2 = and also gives the correct result when p 2 ^ 0. 
F(pi,p 2 ) that satifies these requirements is 




Here, pi,p 2 are 4- vectors whose first three components are the components 
of the momemtum and the fourth component is iE/c, E = (p 2 c 2 + m 2 c 4 )' 5 . 
F{pi,p 2 ) is an invariant and when p 2 = 0, F = 7it>i. The result for F(pi,p 2 ) 
given by Eq.(3) can also be written as 

F(pi, P 2) = 7i7 2 c[(/fi-/f 2 ) 2 -(/fiX/f 2 ) 2 ]- 5 (4) 

Here, /3i,/3 2 are vectors in 3-space corresponding to the velocities of the par- 
ticles divide by c. 

3 The f(x,p) distribution and the scattering 
rate SN 

Let us now treat the case where the particles are contained within a bunch 
and their distibution is given by f(x,p) where Nf(x,p) is the number of 
particles in d 3 xd 3 p. N is the number of particles in the bunch, all particles 
have the same rest mass m and 

J d 3 xd 3 p f(x,p) = 1 

Let 5N be the number of particles with momentum pi in d 3 p± and space 
coordinate x in d 3 x which are scattered by the particles with momentum p 2 
in d 3 p 2 which are also in d 3 x, in the time interval dt , into the solid angle 
dfl' corresponding to the direction p[. Then SN can be obtained using the 
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same procedure used in obtaining Eq.(2), provided one knows that d 3 p/j and 
f(x,p) are invariants, which is shown in section 5. 5N is given by 

5N = N 2 adn' d3pi ^ f( x , pi)f(x, p 2 )F(p u p 2 )d 3 xdt 

7i 72 

F{Pi,P2) = 2 ( 5 ) 

One may note that the right hand side of this expression for 5N is an in- 
variant. We will be putting c = 1 except when something may be gained by 
showing c explicitly. 



4 Growth rates for < piPj > 

Growth rates will be given for < PiPj >. where the <> indicate an average 
over all the particles in the bunch. From these one can compute the growth 
rates for the emittances, < >. The advantage due to computing growth 
rates for < p^pj > stems from the observation that if Pi, pj are the components 
of the momentum 4-vector, then p^pj is a tensor in 4-space and so is 5 < 
PiPj >, as will be seen below, where S < PiPj > is the change in < piPj > 
in a time interval dt. The transfornation properties of a tensor can then be 
used to facilitate the transfer of results between two CS. 

In a scattering event, where a particle with momentum pi scatters off a 
particle with momentum p 2 , the momenta will change to p[ and p' 2 . Let Spu 
represent the change in pu in the collision, and similarly for S(pupij). Then 

Spa = Pu - Pu 
SipuPij) = p'uP'ij - PuPij (6) 

Using the scattering rate given by Eq.(5), one can now compute 5 < 

PiPj > 

<S(p liPlj )> = N[ d*x^^f(x,p 1 )f(x,p 2 )F(p 1 ,p 2 ) 

J 7i 72 

adQ'ip'^j - PuPij)dt 
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One may note that 

< Sipupij) >= 5 < puptj > 

and that 5 < puPij > is a tensor in 4-space because of the transformation 
properties given above for the quantities appearing on the right hand side of 
Eq.(7). Eq.(7) is our general result for the growth rates , holds in all CS, 
and can be used for any particle distribution, f(x,p). 

This result can be further simplified by first considering the integral, for 
a given pi,p 2 , 

Cij = J od^'ip'iiP'ij-PiiPij) (8) 

Cij has the transformation properties of a tensor in 4-space as adQ' is an 
invariant. For a given p±,P2, can be evaluated in the CMS ( the center of 
mass CS ) and if the result can be written in terms of 4-vectors and tensors 
in 4-space, then the result in this form, will hold in all CS. The calculation 
of can be simplified by noting that because of the symmetry in p\ and p 2 
we have 

<S(p li p lj )> = <5(p 2 iP2j)> (9) 
and we can define as 

= J (Tdn'-[5(p li p lj ) + 5(p 2i p 2 j)} (10) 
and Eq.(7) can be written as 

<S( Pli Pij)> = N[ d 3 x d3pid P2 f(x : p 1 )f(x : p 2 )F(p 1 ,p 2 )C ij dt(U) 

J 7i 72 

We will now further evaluate by first evaluating for some particular 
values of p±,p 2 in the CMS corresponding to p±,p 2 and and then using the 
tensor properties of to find a result that holds in any other CS. We are 
particularly interested in finding a result in the Rest CS, which is the CS 
which moves along with the bunch. In the CMS, 

P2 = -Pi 

A = -(fi-P2)=Pi 

ql = p [ - pi 

to = ~<h 
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Using q[ = p[ — p[ and q 2 = — <Ti, one can show that 

\{KPuP\j) + $(P2iP2j) = quqij + Aiqij + Ajqu 



(12) 



In the CMS, we introduce a polar coordinate system 9, where 9 is mea- 
sured relative to the direction of pi or A and we assume that a (9, 0) is a 
fumction of 9 only, we can then write 



pi = (0,0, 1)|A| 

(sin 9 cos 0, sin 9 sin 0, cos#)| A| 
q{ = (sin 9 cos 0, sin# sin 0, cos# — 1)| A| 



P'l 



Considering p\,p 2 to be 4-vectors, and A = .5(p\ — p 2 ), qi = p[ —pi, then A, 
qi are also 4-vectors and in the CMS, A 4 = 0, qu = 0. 

Using Eqs(10) and (12), one now finds for in the CMS 



f'TT 

7T / 

JO 



d9a sin 3 9 |A| 2 



( 1 ^ 
1 
0-20 
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To find in the Rest CS or in the Lab CS, we will try to find an expression 
for in terms of the 4-vectors pu,p2i which gives the above result for C, 
in the CMS. The expression that does this is given by 



c^ 



A,, 



7T f W d9a sin 3 9 A 2 [^- - 3 
J o 



AiA,- 



A 2 



+ 



WjWj 

w 2 - 



hJ 



1,4 



]{PU ~P2i) 



Wi = p U + p 2i 



(13) 



where a is the cross section in the CMS. Let us now verify that this expression 
gives the correct result for Cy in the CMS. In the CMS, 

A 2 = |A| 2 + A 2 = |A| 2 

Wi = i = 1, 3 

A; = i = 1,2 A 3 = |A| 



so that Eq.(13) does give the correct result in the CMS. 
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An important further simplification results from the fact that the particle 
motion is non-relativistic in the CMS and also in the Rest CS which moves 
along with the bunch. For RHIC parameters, for 7 = 100, one finds that 
p ~ le — 3 mc. One can then drop the WiWj/W 2 term. Also A 2 = |A| 2 in 
the CMS and in the Rest CS and one can evaluate F(pi,p 2 ) using Eq.(4) as 
F(puP2) — 2c/3 where f3c is the velocity of either particle in the CMS . In 
the Rest CS, one can write 

dj = 7T r d9asm 3 e [|A| 2 % - 3A;Aj] i,j = 1,3 

J 

Aj = ^(Pli-P2i) 

<$(puPij)> = N J d 3 xd 3 p 1 d 3 p 2 f(x,p 1 )f(x,p 2 )2(3cC ij dt 

fie = \A\/m (14) 

Eq.(14) would be a good starting point for computing growth rates for a 
particle distribution more complicated than a simple gaussian. For the case 
of the Coulomb cross section, one can write CV,- as 

C tj = 2n(r /2p 2 ) 2 \n(l + (2p 2 b max /r ) 2 ) [\K\% - SA.Aj] i,J = 1,3 
a = (r /2(3 2 ) 2 /(l-cos9) 2 

r = Z 2 e 2 /mc 2 (15) 
b max is the largest allowed impact parameters in the CMS. 



5 Invariants d 3 p/j and f(x,p) 

In order to eatablish Eq.(5), one needs to know that d 3 p/'~f and f(x,p) are 
invariants. Consider a CS moving with the velocity v with respect to the 
Laboratory CS. Let the coordinates be x,p in the Laboratory CS and x,p in 
the new CS. p,p are related by 



Ps 


= 7o(p s - v E) 


Px 


= Px 


Py 


= Py 


E 


= \jp 2 + m 2 


7o 
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It then follows that 



dfs = lo(dp s - v dE) 

dE = {p s /E)dp s p x ,p y constant 

dp s = 7o(l - v p s /E)dp s 



E = \Jp 2 + m 2 

E = -f (E-v p s ) 
lQ (l-v Q p s /E) = E/E 
dp s dp s 

EE 1 ' 

Thus dps/'-f is invariant under this momentum transformation and also d 3 p/ r ) 
is invariant. 

Now let us show that f(x,p) is an invariant. Since f(x,p)d 3 xd 3 p is an 
invariant, as it gives the number of particles in d 3 xd 3 p, we need to show that 
d 3 xd 3 p is an invariant. Consider the point x,p in some CS. In the moving 
CS where p s = 0, which is moving with the velocity v = p s /E with respect 
to the first CS, 

dp s = l{dp s - vdE) 
dE = {p s /E)dp s p x ,p y constant 
dp s = dps/'j 
7 = 1/Vl - v 2 (18) 

Since dp s /^ = dp s holds for any CS, dp s /^ is an invariant and d 3 p/^ is an 
invariant. One also has r yd 3 x is invariant because of the Lorentz-Fitsgerald 
contraction. We have then that that d 3 xd 3 p is an invariant. 
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